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We investigate the effect of the dipole-dipole interaction on the quantum jump statistics of three 
atoms. This is done for three- level systems in a V configuration and in what may be called a D 
configuration. The transition rates between the four different intensity periods are calculated in 
closed form. Cooperative effects are shown to increase by a factor of two compared to two of either 
three-level systems. This results in transition rates that are, for distances of about one wavelength 
of the strong transition, up to 100 % higher than for independent systems. In addition the double 
and triple jump rates are calculated from the transition rates. In this case cooperative effects of 
up to 170 % for distances of about one wavelength and still up to 15 % around 10 wavelengths are 
found. Nevertheless, for the parameters of an experiment with Hg + ions the effects are negligible, 
in agreement with the experimental data. For three Ba + ions this seems to indicate that the large 
cooperative effects observed experimentally cannot be explained by the dipole-dipole interaction. 



PACS numbers: 42.50.Ct, 42.50.Ar, 42.50.Fx 

I. INTRODUCTION 

Cooperative effects due to the dipole-dipole interaction 
between atoms are of great importance in many fields, 
most recently in the study of possible quantum comput- 
ers based on trapped ions or atoms, and therefore they 
have attracted considerable interest in the literature [lj. 
A sensitive test for such cooperative effects can be pro- 
vided by atoms showing macroscopic light and dark pe- 
riods in their fluorescence. These can occur in a mul- 
tilevel system if the electron is essentially shelved in a 
metastable state, thereby causing the photon emission to 
cease Q. Two or three such systems accordingly show 
three or four periods of different intensity, namely one 
dark period and bright periods with once, twice, or three 
times the intensity of a single system's bright period. The 
dipole-dipole interaction may alter the statistics of these 
periods. In an as yet unexplained experiment with two 
and three Ba + ions 0, 3 a large number of double and 
triple jumps, i.e., jumps by two or three intensity steps 
within a short resolution time, was observed, exceeding 
by far the value expected for independent atoms. The 
quantitative explanation of such large cooperative effects 
for distances of the order of 10 wavelengths of the strong 
transition has been found difficult 

siania ex- 

periments with other ions showed no observable coop- 
erative effects El i m particular none were seen for 
Hg+ for a distance of about 15 wave lengths p|. More 
recently effects similar to Ref. were found in an ex- 
periment with Ca + ions [Tij l in contrast to a comparable 
experiment [151 ]. A different method for observing the 
dipole-dipole interaction of two V systems was proposed 
in Ref. O. 

The effect of the dipole-dipole interaction for two V 
systems was investigated numerically in Ref. |17| and an- 
alytically in Ref. [18| and shown to be up to 30 % in 
the double jump rate compared to independent systems. 
However, the systems used in the experimental setups of 
Refs. [H Eli E| cannot be described by a V system so that 



a direct comparison between theory and experiment was 
not possible. For this reason the present authors have in- 
vestigated cooperative effects for two other systems |2(1 . 
namely a D shaped system modeling the Hg + ions used 
in Ref. 01 an d a four-level system modeling the Ba + 
ions of Refs. 0, Q . For two D systems cooperative effects 
in the same order of magnitude as for the V systems were 
found for ion distances of a few wavelengths of the laser- 
driven transition. For larger distances practically no ef- 
fects where found, in agreement with the experiments 
[l3| and with the results of Ref. [2l|. In contrast, only 
negligible effects for arbitrary ion-distances were found 
for two of the four level- systems. Although this result 
contradicts the findings of Refs. |||| a direct quantita- 
tive comparison with the experiments was not possible 
since explicit experimental data were only provided for 
three Ba + ions. 

The aim of this paper is to narrow this gap by investi- 
gating three dipole-interacting three-level systems in a V 
configuration and in a D configuration (see Figs. ^ and 
|2J), respectively, and to compare the results with those for 
two such systems. For three system this becomes much 
more complicated since one has to deal with 729 x 729 
matrices, and in order to do this we use group theoretical 
methods to exploit the symmetry of the problem. 

We calculate the transition rates between the different 
intensity periods for both systems. Cooperative effects 
are found to increase by a factor of 2 in the first order 
terms in the interaction parameter C3 when compared 
to two of either systems. This results in transition rates 
up to 100 % higher than the rates for independent sys- 
tems. We also calculate the double and triple jump rates 
for both systems. Here the cooperative effects are even 
larger. 

A full description of the Ba + experiment 0, would 
require the treatment of three of the four-level systems 
of Ref. |2fJ. However, here we will restrict ourselves to 
the three-level systems, since this reduces the complex- 
ity of the calculation considerably. Also, the similarities 
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FIG. 1: Three-level system in V configuration. 



between the results for the D system and the four-level 
system pointed out in Ref. [2(| seem to allow to draw 
conclusions on the cooperative behavior of three four- 
level systems from the results presented here. Namely, 
the increase of cooperative effects is not strong enough 
to yield significant effects for three four-level systems. 

Section^deals with the main assumptions of the mod- 
els. In Section Hill the methods for the calculation of the 
transition rates first for the V systems and afterwards 
for the D system are explained. In Section Hvl the results 
of the calculations are presented, namely the transition 
rates between the different intensity periods. Finally in 
Section[V]the double and triple jump rates are calculated 
from the transition rates. The results are discussed and 
compared with those of two three-level systems. 



II. DIPOLE-INTERACTING THREE-LEVEL 
SYSTEMS 

In the following we investigate three dipole-interacting 
three-level systems both in a V-type and in a D-type con- 
figuration as shown in Figs.nand|21 For the V system the 
Rabi frequencies 1^2 and fl 3 and the Einstein coefficient 
A 3 satisfy 



n 3 ,A 3 > ft 2 



(1) 



so that the single system can show macroscopic light and 
dark periods. The D system exhibits the same property 
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FIG. 2: Three-level system in D configuration with fast tran- 
sitions (solid lines) and slow transitions (dashed lines). 



if the condition 



Q 3 ,A 3 ->a x ,a 2 



(2) 



for the Einstein coefficients and the Rabi frequency is 
fulfilled. We assume the three atoms to be at fixed posi- 
tions forming an equilateral triangle, in agreement with 
the experimental setups. Furthermore, for simplicity, the 
direction of the laser beams are assumed to be perpen- 
dicular to the plane of this triangle. 

The Bloch equation can be written in the form p2^ | 



P : 



H cond p - pHl 



cond 



n{ P ), 



(3) 



where the conditional Hamiltonian H cond and the reset 
operation IZ(p) for a general three-level system are given 
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(1-3 cos 2 9 M ) 



is the coupling parameter which describes the dipolc- 
dipole interaction between atom k and atom I for the 
transition connected with the Einstein coefficient Aj, 
with Oki being the angle between the dipole moments 
and the line connecting the atoms. The dimensionless 
parameter du) = 2irrki/\j is given by the interatomic 
distance ru multiplied by the wave number 2tt/Xj of this 
transition. The detunings of the lasers are taken as zero. 
By setting either A 1 = A 2 = = C$ = or CL 2 = 
in Eqs. J3J and 101 the Hamiltonians and reset states for 
the V systems and the D systems, respectively, are ob- 
tained. For simplicity it would be preferable to have the 
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FIG. 3: Geometry of the atoms in the trap. The arrows 
symbolize the dipole moments. In the picture the angles have 
the values #12 = 7r/2, 6*23 = 7r/6, and #31 = 5tv/6 leading to 
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same coupling parameters for each pair of atoms (i.e., 
C$ = Cj). This would be the case if the angle be- 
tween the dipole moments and the line connecting two 
atoms were the same for all pairs of atoms. However, 
the arrangement of the atoms in the trap makes this im- 
possible, as is illustrated in Fig. [3] The atoms form an 
equilateral triangle (i.e., ru = r) with the laser beams 
perpendicular to the plane of this triangle and the dipole 
moments aligned by a magnetic field in a direction in this 
plane. In this situation, the same value of the coupling 
constants can only be achieved for two of the three pos- 
sible pairs of atoms. However, in spite of this we will 
assume Cu? = Cj because this case leads to maximal co- 
operative effects and can be seen as a limiting case for all 
other possible configurations. The reset state can then 
be written as a sum of density matrices of pure states 



^(p) = ^{(A J +2ReC J ) J R 1 j) P i? i j)t 



therefore use a basis that is adapted to the symmetry 
group S3 of permutations of three particles. On the sub- 
space spanned by the product states with all three atoms 
in different states the irreducible representations of the S3 
are the two one-dimensional representations mentioned 
above and another two equivalent two-dimensional rep- 
resentations. This leads to the states 

| Sl23 ) = -^(|l)|2)|3) + |2)|3)|l) + |3)|l}|2) 



|1)|3)|2) + |2)|1)|3) + |3)|2)|1)), (10a) 
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in the case where all three atoms are in different states. 
For the remaining states one then easily gets for i,j = 
1,2,3, i^j, 
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with 
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In the case of two systems, it was convenient to use a 
Dicke basis, i.e., a basis consisting of the symmetric and 
antisymmetric linear combinations of the product states. 
Generally speaking, this means using a basis which is 
adapted with respect to the symmetry group S2 of per- 
mutations of two atoms. The symmetric and antisym- 
metric states correspond to the irreducible representa- 
tions of this group. For three three-level systems, we 



if two atoms are in the same state and 

\g) = |1)|1)|1), |e 2 ) = |2)|2)|2), |e 3 ) = |3)|3)|3) (12) 
if all three atoms are in the same state. 

III. TRANSITION RATES 

For the calculation of the transition rates, we carry 
over the methods that have already been used for the 
description of two di pole - interacting V systems and D 
systems, respectively fl8t l20| . 

For both types of systems, the configuration decouples 
into four independent subspaces if one neglects the small 
parameters (i.e., VI2 = for the V systems and A\ — 
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A 2 = for the D systems) 

5 = {\e 2 )}, (13a) 

51 = {|S122>, |&122>) |C122>, |»322>, l&322>, |c 3 2 2 )}, (13b) 

5 2 = {|s 2 n>, I&211), |c 2 ii), IS123), |ai 23 ), I&123), 

I cm), \d123), \e12a), \s233), \km), \c233)}, (13c) 

53 = {\g), \s311), \hii), \cau), 

\s133), I&133), |ci33, |e 3 > } (13d) 



where p SSt i is the quasisteady state in subsystem Si. As 
a Liouvillean of Bloch equations, Co has an eigenvalue 
corresponding to the quasisteady states. The other eigen- 
values have negative real parts of the order of and A3. 
While £ip SSj i is a superposition of just the eigenstates for 
nonzero eigenvalues of Co in the case of three V systems 
this is not true for three D systems, which makes it nec- 
essary to discuss the two cases separately. 



in analogy to the case of two of either systems. In a 
period of intensity 7$ , the density matrix of the system is 
mostly in subspace Si [2{| . The transition rates will thus 
be calculated by using a density matrix in one particular 
subspace and then the rate of build-up of population in 
another subspace will be determined. 

Taking a state po.i in one of the subspaces <Sj at a time 
to we calculate the state after a time to + At in pertur- 
bation theory with respect to the small parameters. The 
time interval used here should be long in comparison to 
the mean time between the emission of two photons but 
short in comparison to the length of the intensity periods, 



Ar^nr 1 < At < n? 1 



A3 \ fig 1 < At<Af 1 ,A- 1 



(V system), 
,, (D system) . (14) 



For the calculation the Bloch equation is written in a 
Liouvillean form 



A. Three V systems 



For the V systems, C\p ss ^ consists only of coher- 
ences between the subspace Si and the neighbouring sub- 
spaces, since C\ describes the coupling due to the weak 
laser (with Rabi frequency O2) in this case. The zero- 
eigenvalue subspace of Co, on the other hand, is spanned 
by the quasisteady states p ss ,i- Therefore, C\p ss ^ has no 
components in the zero eigenvalue subspace of Co in the 
case of V systems. The other eigenvalues all have nega- 
tive real parts of the order of A3 and fl3 . Therefore the 
integrand in Eq. i|16[l is rapidly damped which allows us 
to extend the upper integration limit to infinity. This 
yields 



p = Cp={C (A 3 ,C 3 ,n 3 )+C 1 }p 



(15) 



where C\ serves as the perturbation dependin g on H2 or 
A\, A%, Ci, and C2, respectively. We then get [Tg 



f At c 

p(t + At-p ssA ) = p ssA + / dre ° T Cip 

Jo 



(16) 



p(t Q + At; p SSj4 ) = p sM + (e - £ ) 1 Cip ss ^ , (17) 

independent of At . 

From the Bloch equations J3J we get the exact relations 
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= — f^Im 



2(Sll2|p|si 22 ) - (&112 I/0I&122) - (Cn 2 |p|ci22) - V / 3(si22|/0|e2) + V2(si23|p|s223> (18c) 
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dt S MA 
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Together with Eq. (|17f) this allows us to calculate the 
transition rates as 



Pu 



dt 



P=p(*o+At;p as ,i) 



(19) 



Note that = for \i — j\ > 2 so that no direct, i.e., 
instantaneous, double jumps occur. 



B. Three D systems 

In the case of D systems, C\ describes spontaneous 
emission due to the Einstein coefficients Ai and A 2 . 
Therefore C\p ss ^ consists of density matrix elements 
(xi\p\xj) where both states \xi) and \xj) lie in the same 
subspace Si. It is thus a superposition of eigenstates of 
Co with zero as well as nonzero eigenvalues. We write 



C\p as ,i = ^ OlijPssJ + p, 
3=0 



(20) 



where p contains the contributions from the eigenstates 
for nonzero eigenvalues of Co- The coefficients ctij are 
calculated by means of the dual eigenstates p z ss [2fl |. 

= Tr(/^|£ip ss>i ). (21) 

Inserting Eq. (|20|l into Eq. i|16|) one obtains 

3 

p(t Q + At) = p SfM + CHjfhmj&t + (e - Co^p . (22) 

3=0 

The last term is much smaller than the preceding term 
and can be neglected j2jj. The coefficients can then 
be interpreted as the transition rates between the sub- 
spaces Si and Sj , 



Pij = Oti 



C. Group theory 



(23) 



For the calculation of the transition rates for both V 
systems and D systems it is necessary to calculate the 
quasisteady states p ss ,i, i-e., to solve the linear equation 

CoPss = 0. (24) 
In addition, for V systems the first order term 

= (e- Coy 1 Cip ss , i 
of Eq. (|17|l must be calculated, which was done by solving 



two different properties of Co that make it possible to re- 
strict these equations to smaller subspaces, which reduces 
the calculation effort considerably. First, Co is indepen- 
dent of the small parameters (A\, A 2 or 2 ), which means 
that there is no coupling between the four subspaces of 
Eq. (fH5jl . Thus there exist 16 subspaces TZij, each con- 
sisting of the density matrix elements 



with 



€ Si and \ijj) € S 



(26) 



respectively, which are invariant with respect to Co- In 
addition the conditional Hamiltonian -ff C ond and the re- 
set state 7Z(p) and therefore also Co are invariant under 
the exchange of atoms, as can be seen from Eqs. (0} and 
(0. Hence subspaces which consist of all density matrix 
elements which belong to a particular irreducible rep- 
resentation of £3 are also invariant with respect to Co- 
Since the density matrix elements form a representation 
of 53 which is a tensor product of twice the represen- 
tation spanned by the Dicke basis of Eq. IjlOjl the new 
irreducible representations are easily found. The density 
matrix elements 

|s Q )(s /3 |, \a a ){a^\, 

l(K)(bp\ + \c a )(c p \), l -(\d a ){dp\ + \e a ){e p \), (27) 

\{\b a ){e \-\c a ){d p \), \{\e a ){bp\-\d a ){cp\) 

belong to the symmetric representation, the elements 
\s a ){ap\, \a a ){sp\, 
l(K)(cp\-\c a )(b p \), I(|d Q )( e/3 |-|e Q )(^|), (28) 

\ (\ba)(dp\ + \c a ){ep\) , - (KXM + \e a ){ Cf )\) 

belong to the antisymmetric representation, and the 
remaining 24 possible linear combinations form two- 
dimensional representations. Here a and (3 are one of 
the subscripts of the Dicke states. By transforming the 
Liouvillcan Co into this new basis each of the 16 invari- 
ant subspaces TZi,j is in itself decomposed into three in- 
variant subspaces connected to the elements belonging 
to the symmetric, antisymmetric, and two-dimensional 
representations, respectively. For the calculation of both 
the quasisteady states p ss .i and the transition rates for 
the V systems only the symmetric subspaces are needed. 
For the latter this can be seen from Eq. (|18fl . With these 
two simplifications the dimension of the linear system of 
equations needed for the calculation reduces considerably 
(namely to a maximum of 20 for the calculation of P23 
and p 32 ). 



r C 1 ) r 
CoPi = Cip s 



(25) 



IV. RESULTS 



Eqs. (|24|l and Ij25(l are linear equations for the 729 matrix 
elements of p SSi i and p^\ respectively. Luckily there are 



The transition rates for the V systems can now be cal- 
culated according to Eqs. (|T§|l . l(TB|) . and (jT%|) . The result 
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(29c) 
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(29d) 
(29e) 
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to first order in C3. While for poi and pio this is also 
the exact result to all orders, the higher order terms for 
the other four transitions are too complicated to be given 
here. The zeroth order terms in Eqs. (|29|l are those one 
would expect for independent atoms (namely the rates 
Piq and poi for single V system multiplied by a factor 
1, 2, or 3). For the first order terms it is important to 
note that the single systems interact via C3 only if they 
are in a light period. Therefore the rates poi and pio are 



independent of C3 while P12 and P21 have the same first 
order term as the corresponding rates for two V systems 
(in the intensity period I2 the three V systems behave 
like two V systems in the period I2 plus an additional 
nonintcracting system). In the rates P23 and P32 the first 
order term is just twice the first order term of P21 and 
p±2- This surprising property is due to the simplicity of 
the quasisteady state p S s.3, namely all diagonal elements 
of this state have the same first order dependence. Fig.0] 
shows the transition rate P32 for three V systems to first 
and to second order in C3. The first order rate becomes 
negative for distances of about one-half to three quarters 
of a wavelength of the strong transition. By looking at 
the second order rate one can see that is an artefact of 
the approximation. The rate with the dipole interaction 
included shows deviations of up to 100% from the rate 
for noninteracting atoms for distances of somewhat more 
than a wavelength A3. 

By use of Eqs. J2U and J22Jl the transition rates 

for three dipole interacting D systems were also calcu- 
lated, with the result 



p 01 =3A 1 , p 12 = 2A 1 , P23 = A U 
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P21 



P32 



[Al 



20|)2 



2A 2 0§ 

^§(|C 3 | 2 + 2 J 4 3 ReC , 3 ) ~ A\ + 20| 



A 2 nl(Al 



20|) 



1 - 2Re C; 



At 



3A 2 nj [(A 



202)2 +A 2(|C 3 | 2 



(Al + 201)2 
+ 2A 3 ReC 3 )] 



0(C 2 3 ), 



(A|+20|)[(A1 



202)2 + 3A 2 3 (\C 3 



3A 2 



Al 



2"3 

202 



1 - 4Re C 3 



(At 



A 3 3 



20§)2 _ 



2 + 2A 3 ReC 3 
+ 0{Cj). 



2A2[|C 3 |2|A 3 
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(30d) 



Compared to two D systems the transition rates show 
the same behavior as explained above for the three V 
systems. This is not surprising as the quasisteady states 
are identical and as the D systems also only interact via 
C3 when they are in a light period. Fig.[S]shows the exact 
transition rate P32 compared to the interaction free case. 
For distances of about a wavelength, p 3 2 deviates up to 
75 % from the rate without interaction. The first peak at 
about 0.7 wavelengths even reaches a maximum of seven 
times the rate for independent atoms. For such small 
distances, however, one would have to check the validity 
of the model (namely, that in a particular intensity period 
most of the population is in a specific subspace). Also 
one must keep in mind that all the experiments cited here 
were performed at greater ion distances. 



V. DOUBLE AND TRIPLE JUMP RATE 

The physical quantity investigated in the experiments 
of Refs. |E H HI [H HI is the double jump rate. This 
is the rate at which jumps between periods of intensities 
that differ by twice the intensity of a single system occur 
within a small time interval. In Ref. |l8j the double jump 
rate has been expressed in terms of the transition rates 
Pij for two dipole-interacting V systems. The same will 
be done here for three systems. As one can calculate di- 
rectly from Eqs. (|19fl and l|23|l there are no direct double 
jumps (i.e., py = for \i — j\ > 1). A double jump is 
therefore defined as two successive jumps in the same di- 
rection which occur within a time which is smaller than 
a time window T m so that they cannot be resolved. As 
there are four periods of different intensity in the fluores- 
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FIG. 4: Transition rate P32 for three dipole-interacting V 
systems plotted versus the interatomic distance r in units of 
the wavelength A3 of the strong transition. Solid line: P32 up 
to second order in C3. Dashed line: first order. Dotted line: 
independent systems. Parameter values are A3 = 2- 10 s s -1 , 
Q 3 = 5- 10 7 s _1 , and Q 2 = 10 4 s _1 . 
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FIG. 5: Transition rate P32 for three dipole-interacting D 
systems. Dashed line: independent systems. Parameter 
values are A\ = Is -1 , Ai = Is -1 , A 3 = 2 ■ 10 s s" 1 , and 
fi 3 = 10 7 s _1 . 



cence of three three-level systems, there are also four dif- 
ferent possibilities for double jumps: From intensity zero 
to double intensity, from single intensity to threefold in- 
tensity, and vice versa. Therefore the whole double jump 
rate «dj is the sum of rates for the four different possible 
double jumps, 

_ „20 1 „31 , „13 1 „02 (ot \ 

We first derive the rate for jumps from zero to double 
intensity. Each period of zero intensity ends with one of 
single intensity. The probability that the latter period is 
shorter than T m is given by 

PT 1<Tn , = 1 - e -< Pl0+, ' ls > T » . 

The branching ratio for the following period to be of 
double intensity ispi2/(f>io +Pi2)- With the mean num- 
ber of intensity periods Ii per unit time denoted by m 

0.5 r-p 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 \ 



0.4: 




I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I I 

0123456789 

r[As] 

FIG. 6: Double jump rate noj for three dipole-interacting V 
systems. Solid line: noj up to second order in C3. Dotted 
line: independent systems. Time window T m = 10 -3 s. Other 
parameter values as in Fig. 0] 



the rate n^j is given by 

n 02 = P12 / ! _ e -( P10+Pl2 )T m \ (32) 

DJ P10+P12 V / 
Analogously one finds 

n 3i = n3 P21 ( 1 _ e -( P 2 1+P23 )T m \ _ (33) 
P21+ P23 V > 

The remaining two rates are a little bit more compli- 
cated as the periods of intensity 1\ and I2 can be followed 
by a period with either higher or lower intensity. The 
rates nj^j and n^j have thus to be supplemented with 
the branching ratios p2i/(P2i +P23) and pi 2 /{pio +P12), 
respectively, yielding 




0123456789 

r[As] 

FIG. 7: Double jump rate noj for three dipole-interacting 
D systems. Dashed line: independent systems. Time window 
Tm = 5 • 10 -3 s. Other parameter values as in Fig.0 
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Using the the relations 
Pw 



n 



Pw +P12 



-m, n 3 



P23 



P21 + P23 



-ri2 



and 



Tl2 = 



P12 



-n\ + n 3 , n\ = n + 



P21 



P1Q+P12 P21+P23 
the double jump rates can be simplified to 



(36a) 



-712 (36b) 



02 



20 
'DJ 



ni 



P10P12 



and 



'D.l 



n 



D.T 

m 



(pw +P12) 



P12P23 



l-e 



-(P10+P12)T„ 



(37) 



(38) 



(P21 +P23)(PW + Pl2) 



-(P21+P23)T m 



We denote the mean durations of the intensity periods 
by Ti and note that 



T = — , T 1 = 1 , T 2 = 

P01 P10 + P12 



1 



T 3 = 



P21 + P23 



In addition they fulfill 



1 

P32 ' 

(39) 



(40) 



i=0 



The averaging window T m is much smaller than the mean 
durations of the intensity periods. Therefore the expo- 
nential can be expanded and with Eq. lpH|) one gets 



Puipw +P23) 



PlO + Pl2 

Using Eqs. lj3li]). l|39")l. and ijHU we finally obtain 
PoiP2iP32(poi +P12) 



(41) 



P2lP32(P01 +Pw) +P01P12(P23 + P32) 



T m (42) 



as the double jump rate for three of either three-level 
systems. A similar calculation yields for the triple jump 
rate 



"-TJ 



P01PWP12P21P23P32 



P2lP32{Pm + Pw) + P0lPl2(P23 + P32 



(43) 



Note that the defining time window T m enters quadrat- 
ically in this case. Figs. El and show plots of tidj for 
the V systems and the D systems, respectively, whereas 
Figs. |H1 and show plots of the triple jump rate titj for 
both systems. For the D systems the exact values for the 
Pij are used whereas for the V systems only the expanded 
expressions up to second order in C3 are used since P23 
and P23 could not be calculated exactly for the V sys- 
tems. For the V systems there are cooperative effects of 
up to 110% for the double jump rate udj and 170% for 
the triple jump rate nxj for distances of somewhat more 




FIG. 8: Triple jump rate titj for three dipole-interacting V 
systems. Solid line: jitj up to second order in C3. Dotted 
line: independent systems. Parameter values as in Fig. |S| 
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FIG. 9: Triple jump rate nxj for three dipole-interacting 
D systems. Dashed line: independent systems. Parameter 
values as in Fig. [T] 



than a wavelength of the strong transition. For the same 
distance range the D system shows cooperative effects of 
up to 150% for both 71dj and tjtj- The first peak at 
three quarters of a wavelength reaches 16 times the value 
for independent systems for both rates. For distances of 
about 10 wavelengths cooperative effects of 15 % are still 
present for both systems. In the case of the D system, 
which models the level confi gur ation of the Hg + ions used 
in the experiments of Refs. |13ll26T |. large cooperative ef- 
fects only appear if the Rabi frequency is smaller than 
the Einstein coefficient A3. So, for the experimental pa- 
rameters (i.e., O3 > A3 and r/X^ w 15) the effects are 
negligible, in agreement with the experimental results. 



VI. CONCLUSIONS 

We have investigated the effect of the dipole-dipole 
interaction on three three-level systems showing macro- 
scopic light and dark periods in their fluorescence. This 
was done for the V and the D configuration, respectively. 
The latter models the effective level configuration of the 
Hg + ions in the experiments of Refs. |13j,|26j. We have ex- 
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plicitly calculated the transition rates between the differ- 
ent intensity periods for both configurations. In addition, 
the double and triple jump rates have been derived from 
these transition rates. Both systems show the same first 
order dependency on the coupling parameter C3, leading 
to an enhancement in the cooperative effects by a factor 
of 2 for the transition rate p^2- This leads to cooperative 
effects of about 100% compared to the value for inde- 
pendent systems for interatomic distances of somewhat 
more than a wavelength of the strong transition. For 
the double and triple jump rates even larger cooperative 
effects can be seen. For three D systems the first peak 
at about three quarters of a wavelength is seven times 
higher for P32 and 16 times higher for tidj and riTj than 



for independent atoms. 

Although we did not treat the four-level system of 
Ref. [2(j here, which models the Ba + ions of Refs. [H,Qt, 
it is still possible to arrive at some conclusions on this 
experiment from our results for the three level systems. 
As was pointed out in Ref. [2fij . the results for two D 
systems and two four-level systems are very similar, in 
particular in their first order term in C3. It is therefore 
very likely that the cooperative effects for three four-level 
systems are also only enhanced by a factor of about 2, 
and since the effects for two four-level systems were al- 
ready negligibly small one can expect a similar behavior 
also for three of such systems. 
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